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ABSTRACT  (JUAlm  WSPEOm,  3 

General  models  describing  the  interactions  between  a  pair  of  piezoceramic  patches  and 
elastic  snhstrnctures  consisting  of  a  cylindrical  shell,  plate  and  beam  are  prc'sented.  In  each 
case,  the  manner  in  which  the  patch  loads  enter  both  the  stro?)g  and  weak  forms  of  the  time- 
dependent  structural  eciuations  of  motion  is  described.  Through  force  and  moment  balancing, 
these  loads  are  then  determined  in  terms  of  material  properties  of  the  patch  and  substructure 
(thickness,  (dastic  jiroperties,  Poisson  ratios),  the  geometry  of  the  patch  placement,  and  the 
voltages  into  the  patches.  In  the  case  of  the  shell,  the  coupling  between  bending  and  in¬ 
plane  deformations,  which  is  due  to  the  curvature,  is  retained.  These  models  are  sufficiently 
general  to  allow  for  jjotentially  different  patch  voltages  which  implies  that  they  can  be  suit¬ 
ably  employed  when  using  piezoceramic  patches  for  controlling  system  dynamics  when  both 
extensional  and  bending  vibrations  are  [)resent. 
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1  Introduction 


The  use  of  piezocerainic  elements  as  sensors  and  actuators  lias  burgeoned  in  the  last  several 
years  in  applications  ranging  from  the  measurement  and  damping  of  vibrations  in  large  flexi¬ 
ble  structures  to  the  control  of  noise  in  structural  acoustics  settings.  Their  utility  as  sensors 
derives  from  the  property  that  when  the  element  is  subjected  to  a  mechanical  strain,  a  voltage 
proportional  to  the  strain  is  produced.  Conversely,  they  also  exhibit  the  phenomenon  that 
an  applied  polarization  voltage  across  the  unconstrained  element  produces  in-plane  mechan¬ 
ical  strains  in  the  material.  Because  of  these  properties,  piezoceramic  elements  have  found 
increasing  success  both  as  sensors  such  as  strain  gauges  and  accelerometers  and  as  distributed 
actuators.  Their  success  as  actuators  is  augmented  by  the  fact  that  they  can  be  used  to  di¬ 
rectly  control  local  vibrations  without  applying  rigid  body  forces  and  torcpies,  and  due  to  their 
distributed  nature,  they  are  less  prone  to  spillover  effects  in  many  control  strategies.  More¬ 
over,  the  piezoceramic  elements  or  patches  are  inexpensive,  lightweight,  space  efficient  and 
can  be  easily  shaped  or  bonded  to  a  variety  of  surfaces.  Hence  a  large  mmiber  of  patches  can 
be  used  to  sense  and  control  without  significantly  changing  the  mass  or  dynamic  properties 
of  the  system. 

In  order  to  obtain  optimal  results  with  the  piezoceramic  elements  or  patches  in  sensing 
and  control  applications,  it  is  necessary  to  have  accurate  models  of  the  nuM  lianics  of  induced 
strain  actuation.  This  modeling  also  provides  knowledge  of  tlu'  physical  limitations  of  tlu' 
piezoceramic  patches  as  actiiators  in  various  settings.  Detailed  modt'ls  have  l)e('n  d(>veK)ped 
for  piezoceramic  patch  interaction  with  Euler- Bernoulli  beams  [2.  2.  1,  5.  10]  and  thin  i)lates 
[6,  12].  Because  many  of  the  initial  applications  of  piezoceramic  ehmuMits  were  in  settings 
involving  the  sensing  and  control  of  bending  deformations  (these  vibrations  arc'  dominant  in 
many  low  frequency  vibration  and  noise  control  problems),  most  of  these  models  concc'iitrate 
on  patch  configurations  which  excite  pure  bending  motion  of  the  substructure  with  more 
limited  discussions  of  pure  extensional  excitation.  It  was  not  until  [10]  that  a  model  was 
developed  which  provided  for  simultaneous  excitation  of  both  bending  and  extensional  de¬ 
formations  in  an  Euler-Bernoulli  beam.  One  motivation  for  developing  such  a  model  is  the 
obsc'rvation  that  in  complex  coupled  systems,  in-plane  vibrations  with  small  dis])lac('ments  can 
have  large  in-jdane  energy  levels  due  to  the  property  that  beams  are  much  stiffer  in  extension 
than  in  bending.  This  in-plane  energy  can  then  couple  into  flexural  vibrations  at  structural 
discontinuities  such  as  joints,  thus  necessitating  the  control  of  both  bending  and  extensional 
vibrations  in  such  structures.  As  determined  by  Fuller  et  al  [8]  through  experimental  work, 
simultaneous  reductions  in  both  flexural  and  extensional  deformations  in  a  beam  can  be  ob¬ 
tained  through  the  use  of  asymmetric  pairs  of  piezoceramic  actuators  and  sensors  in  adaptiv(' 
control  schemes,  and  the  analytic  work  in  [10]  was  a  first  step  toward  developing  a  moilel 
which  could  be  used  in  further  such  control  settings.  In  that  work,  force  and  moment  balanc¬ 
ing  were  used  to  determine  ex|>ressions  for  the  moments  and  strains  induced  by  tlu'  activation 
of  a  single  piezoceramic  patch  which  was  bonded  to  an  I’hih'r-lh'rnoulli  beam. 
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In  addition  to  lieains  and  plates,  thin  clastic  shells  an'  often  used  to  tlescrilx'  various 
strncttiral  components  as  well  as  when  modeling  the  coupling  hetween  structural  vihrations 
and  their  radiating  or  receiving  acoustic  fields.  For  example,  the  t  ransmission  of  sonnd  through 
an  airplane  fuselage  (hie  to  low  fre(]nency,  high  amplitiuh'  I'xtc'rior  acoustic  ii<'lds  can  Ik* 
modeh'd  by  a  vilirating  thin  cylindrical  shell  which  is  coniih’d  to  an  interior  acoustic  pressure 
fi(’ld  [f)].  In  order  to  optimally  control  the  inlc'rior  nois«'  via  piezoc('ramic  patch  actuation, 
one  first  needs  to  accurately  nujdel  tin'  interactions  hetwc'c'ii  the  patclu's  and  the  shell.  This 
raises  modeling  issues  which  differ  frcim  those  encoimlerc'd  in  tin-  beam  and  plate  analyses  in 
that  the  in-plane  and  bending  vibrations  arc'  coupled  in  the  cylindrical  she'll  du('  to  curvatuix' 
effects. 

Analytical  models  describing  piezoc'eramic  patch/c\  lindrical  sht'll  inti'ract  ions  have'  primar¬ 
ily  ht'en  based  on  layere'd  shell  tlu'ory  [11.  1!)]  or  the'  use'  of  Hat  plate'  pie'zoce'iamic  coujding 
re'sults  when  determining  the  re'sulting  loading  on  the  she'll  [lol  bi  the'  first  case'  it  is  assunu'd 
that  the  pie'zoceramic  material  make's  np  an  entire'  laye'r  of  the  e'lastic  structure  and  hence 
this  model  is  of  limite'd  use  wdien  considering  small  patche's  as  actuators.  When  using  the' 
flat  plate  the'ory.  it  is  assumed  that  the  patch  dinu'nsions  are'  small  in  comparison  with  the 
cylinder  radius,  (’urvature  propertic's  are  then  ne'glected  whe'ii  moeh'ling  the'  coupling  be'twe'en 
the  patch  and  shell  and  de'termining  the  loading  due  to  activation  of  the  jeatch. 

In  this  work,  we  present  general  models  for  the  interactions  be'twe'e'ii  a  pair  of  piezoceramic 
patches  and  e'lastic  sid)strnctures  consisting  of  a  be'am.  plate  or  thin  she'll.  In  the'  case  of  a 
shell,  the  jjatclu's  are  assumed  to  be  curved  and  the  coei))ling  betwe'en  bending  and  in-plane 
deformations,  which  is  due  to  the  curvature,  is  retaiiu'd.  The  teehnie|ues  used  to  deve-loj)  the 
shell/patch  interaction  model  are  also  used  to  develop  general  models  describing  the  moments 
and  forces  which  are  generate'd  by  the  activation  of  piezoceramic  patches  which  have  here'll 
boneh'd  to  a  flat  plate  or  beam.  These  models  dilfer  from  those  in  [h.  12]  and  [2.  3.  4.  o]  in 
that  the'v  allow  for  different  voltages  into  the  individual  patche's  thus  admitting  the'  analysis 
of  simultane'ous  excitations  of  both  bending  and  extensional  comj)onents.  In  the  case  of  the' 
be'am.  the  motlel  is  slightly  more  general  than  that  in  [10]  sinee'  it  is  de'rive'd  for  twee  active 
patche's.  Me'iice  it  is  appropriate  for  moeh'ls  in  which  both  patches  actuate  with  potentially 
dilfe'ient  voltages  into  the  patclu's.  In  the  case  of  one  actuating  patch,  howeve'r.  tlie  model 
reduce's  to  that  in  [10].  From  a  control  perspective,  the'se  moeh'ls  are  important  since  the'v 
provide  for  greater  latitude  in  eh'signing  control  slrate'gie's  involving  the  eise'  of  pie'zoceramic 
e'lements  to  affect  both  the  bending  and  extensional  propertie's  of  a  structure. 

As  a  prelude'  to  the  development  of  the  patch  inte'raction  mode'ls,  e'(|uations  of  motion  for 
the  uneh'rlying  substructure's  are  pre.sente'd  with  special  atte'ution  j)aid  to  the  contributions 
due'  to  ('xternally  apjdie'd  moments  and  forces  since  this  is  where  the  interactions  between 
the  [jatche's  and  sidestructure  occur.  The  analysis  h'ading  np  to  the  structural  e'epiations  also 
motivate's  many  of  the'  te'chniepu's  which  are  use'd  to  eh'velop  the'  patch  interaction  models. 

To  this  ('iid.  a  synopsis  of  the'  derivation  of  the'  strong  form  of  the  time'-de'ix'ndent  Donnell- 
.\lnshtari  thin  she'll  ('e|uations  from  Ne'wtonian  i)rinci|)l«'s  (force  and  moment  balancing)  is 
pre'sente'd  in  Se'ction  2.  A  complete  treatme'iit  of  this  topic  can  be  found  in  [13,  14,  Ki.  17, 
1<S]  and  the'  incbieh'd  discussion  is  limite'd  to  summarizing  that  mate'rial  which  is  ne'e'de'd  for 
de've'loping  the  shell/patch  interaction  mode'l  as  pre'sente'd  in  the  following  se'ction.  The'  choice 
of  the  Donne'll-Mushtari  mode'l  is  for  ease' of  pre'se'utation  and  as  note'd  at  various  points  in  the' 
discussion,  the'  patch/shell  interaction  moeh'l  can  be  e'asily  exte'nde'd  to  higlu'r  orele'r  moeh'ls 


as  warranted  by  the  physical  situation. 

An  inherent  disadvantage  of  the  strong  form  for  the  equilibrium  equations  when  the  ex¬ 
ternal  loads  are  generated  by  piezoceramic  elements  is  the  resulting  presence  of  the  first  and 
second  derivatives  of  the  Heaviside  function  due  to  the  finite  support  of  the  jjatches.  As  a 
result  of  this  as  well  as  other  identification  and  approximation  issues,  we  then  develop  the 
weak  form  of  the  time-dependent  Donnell-Mushtari  shell  equations.  This  is  done  in  more 
detail  since  this  development  is  less  readily  available  in  the  the  literature.  This  forimdation  is 
advantageous  in  many  approximation  schemes,  admits  the  identification  of  discontinuous  ma¬ 
terial  parameters,  and  eliminates  the  problem  of  differentiating  the  Heaviside  function  since 
the  derivatives  are  transferred  onto  the  test  functions. 

The  second  section  concludes  with  a  .synopsis  of  the  strong  and  weak  forms  of  the  Kirchhoff 
plate  and  Euler-Bernoulli  beam  equations.  As  in  the  shell  discussion,  particular  emphasis  is 
placed  on  the  contributions  of  externally  applied  forces  and  moments  since  this  is  where  the 
coupling  between  the  substructure  and  piezoceramic  patches  occurs. 

The  patch  contributions  to  the  cylindrical  shell  equations  are  developed  in  Section  In 
order  to  determine  the  loading  due  to  patch  moments  and  forces,  it  is  \iseful  to  first  express 
them  in  terms  of  the  normal  strains  and  changes  in  curvature  of  the  middle  surface  of  the 
cylindrical  shell.  To  do  this,  the  stress-strain  relations  in  the  patch  and  shell  are  developed 
which  then  allows  the  moment  and  force  resultants  for  the  patch  to  be  formulated  in  terms  of 
midsurface  shell  properties.  The  unknown  normal  strains  and  midsurface  changes  in  ciirvat  ure, 
and  hence  the  patch  moments  and  forces,  are  then  determined  by  moment  and  force  balancing. 
In  this  manner,  the  loading  due  to  activation  of  the  patches  can  be  expressed  in  terms  of 
material  properties  of  the  patches  and  shell  (thickness,  elastic  properties,  and  Poisson  ratios), 
the  radius  of  curvature  of  the  shell,  and  the  voltage  being  applied  to  the  patches. 

In  Section  4,  the  techniques  of  the  third  section  are  tailored  to  composite  structures  consist¬ 
ing  of  piezoceramic  patches  which  are  bonded  to  plates  and  beams.  The  resulting  plate/patch 
interaction  model  is  shown  to  be  equivalent  to  that  of  [12]  in  the  special  case  when  pure 
bending  motion  is  excited  (the  model  in  [12]  was  derived  by  isolating  the  interface  stress  for 
the  system  and  treating  it  as  the  unknown  to  be  determined).  Due  to  its  generality  however, 
our  model  also  allows  for  more  complex  interactions  involving  both  bending  and  extensioi.cl 
components  since  the  voltages  into  the  individual  patches  can  diifer.  As  discussed  earlie’ .  i  he 
beam/patch  interaction  model  reduces  to  the  model  in  [10]  in  the  case  of  one  actuatin  ;  patch 
but  is  slightly  more  general  in  that  it  also  admits  models  in  which  two  patches  an  used  for 
actuating  with  potentially  different  voltages  into  the  patches.  As  with  the  shells,  this  provides 
structure/patch  interaction  models  which  can  be  used  in  various  structural  and  structural 
acoustics  control  settings. 


2  Underlying  Shell,  Plate  and  Beam  Equations 


I  lirougliDut  this  discussion,  wo  consi<l(>r  a  t  Idii  circular  cylindrical  slu-ll  of  radius  /f.  t  liickncss  h 
and  having  the  axial  (X)ordinate  .r  as  shown  in  Figure  1.  1  he  xariahh'  c  ineasuia's  the  distance 
of  a  |>oinf  on  tlie  shell  from  the  <  orresponding  |)oint  on  tin'  middle  surface  (c  =  0)  along  tin' 
normal  to  tlie  middle  surface. 


Figure  1.  The  (’ylindrical  Thin  Shell. 


Strain-Displacement  Relations 

By  combining  Love's  shell  assumptions  witli  the  strain-displacc'im'iil  (’(|ualions  of  three 
dimensional  elasticity  theory,  oik’  obtains  the  strain  ndations 

( j.  —  t  j-  T  „  ^  j- 

=  (2.1) 

where  ( j.  and  ( o  are  normal  strains  at  an  arbitrary  point  within  the  cylindrical  shell  and  ')j.p  is 
the  slu'ar  strain.  Here  cj-.  £o  and  £j-o  are  the  normal  ainl  shear  strains  in  the  middle  surface  and 
Kj.,  h'o  and  r  are  the  midsurface  changes  in  curvature  and  midsurfact'  twist  (sec'  [11],  page  8). 

.\ot<'  that  within  the  framework  of  infinitesimal  elasticity,  the  eciuations  (2.1)  are  exact 
and  in  the  Byrne- Fluggc'-Lur’ye  shell  theory,  these  represc'iit  th<'  c'xact  forni  of  tin'  kinc'inafic 


equations.  In  tlie  Donnell-Mushtari  theory,  one  neglects  the  underlined  terms  z  j R  with  respect 
to  unity  thus  leaving 

f'j.  =  +  ZK^ 

^0  =  ^9  +  -  tie  10 


-  e,.e  +  -  (l  +  ^ 


In  terms  of  the  axial,  tangential  and  radial  displacements  u,  v  and  re,  respectively,  the 
ex|)ressions  for  the  midsurface  strains  and  changes  in  curvature  for  the  cylindrical  shell  are 
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As  before,  the  underlined  terms  are  retained  in  the  Byrne,  Fh'igge  and  Lur'ye  theory 
and  are  discarded  in  the  Donnell-Mushtari  theory.  We  point  ont  that  the  equations  (2.1) 
and  (2.3)  differ  from  those  arising  in  the  theory  of  flat  plates  both  in  the  presence  of  the 
length  differential  RdO  as  well  as  in  the  retention  of  the  strain  terms  tj-  and  eo  (only  bending 
contributions  are  considered  in  the  corresponding  models  of  the  transverse  vibrations  of  a  flat 
plate). 

Stress-Strain  Relations 

To  determine  the  constitutive  properties  of  the  shell,  it  is  assumed  that  the  shell  material  is 
elastic  and  isotropic.  Hooke’s  law  in  conjunction  with  the  assumption  that  the  t  ransverse  shear 
.stre.sses  cTj.-  and  ag,  as  well  as  the  normal  strain  component  are  small  in  comparison  with 
other  stresses  and  strains  (these  conditions  are  part  of  Love’s  third  and  fourth  assumptions) 
then  yields 

£ 

<^r  -  - - 7((’j-  +  t'Vo) 

1  — 

(re  =  - - T(f0  +  /^fr)  (2.1) 


where  <7^  and  ag  are  normal  stresses  aiul  rr^g  and  ag^.  ar<‘  tangential  shear  stresses.  The 
constants  E  and  i/  are  the  Young’s  modidus  and  Poi.sson  ratio  for  the  shell. 


Force  and  Moment  Resultants 


By  integrating  the  stresses  over  the  face  of  a  fundamental  element,  the  force  resultants  can 
be  expressed  as 


Similarly,  the  moment  resultants  are 


'Me  1  ^  ere 


The  orientation  of  the  various  forces  and  moments  are  shown  in  Figure  2.  We  point  out  that 
the  transverse  shear  stresses  and  cr^;  are  used  when  obtaining  the  force  resultants  Qj.  and 
Qe  even  though  they  are  omitted  in  the  constitutive  relations.  This  is  one  of  the  contradictions 
which  arises  in  the  classical  shell  theory. 


Q, 


Figure  2.  F’orce  and  Moment  Resultants  for  the  Cylindrical  Shell. 


In  the  Donnell-Mushtari  theory,  the  underlined  terms  z/ R  are  neglected  in  comparison  to 
unity,  and  the  integrals  are  determined  accordingly  to  y  .d 
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(2.5) 


Similar  expressions  are  obtained  in  the  higher  order  theories. 


Strong  Form  of  the  Donnell-Mushtari  Shell  Equations 


The  equations  of  the  dynamic  equilibrium  of  the  element  are  obtained  by  balancing  the 
internal  force  and  moment  resultants  as  shown  in  Figure  2  with  any  externally  applied  forces 
and  moments.  Let 

<7  =  (lAx  +  90^0  T  </n*n 

and 

m  =  ihjjr  +  dlglg 

denote  the  surface  forces  and  moments  due  to  an  external  field  which  is  acting  on  the  middle 
surface.  Hence  q  and  m  have  units  of  force  and  moment  per  unit  area,  respectively. 
Considering  equilibrium  of  the  forces  in  the  .r,  0  and  directions  yields 
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respectively.  In  the  Donnell-Mushtari  theory,  the  transverse  shearing  force  Qg  is  considered  to 
be  negligible  in  the  second  equation  of  (2.6)  and  is  subsequently  neglected  when  determining 
the  final  equilibrium  equations.  Similarly,  with  0  as  a  reference  origin,  the  balancing  of 
moments  with  respect  to  0,  x  and  c  yields 


dM,  dMgr 

R~r. - h  — - RQt  +  Rule  —  0  , 


dx 
dMg 

~W 


de 

+  —  hfQe  +  Rdir  =  0  . 

dx 

Nre-Ng,-~=0. 


(2.7) 


R 


respectively.  By  referring  to  the  integral  definitions  of  aiid  Moj  -  it  can  be  .seen  tliat 

the  third  e.xpression  in  (2.7)  is  identically  satisfied  due  to  the  syninietry  of  the  stri'ss  tensor. 

Time  enters  the  equilibrium  erpiations  through  the  inert  ial  terms;  hence  for  time-dependent 
problems,  the  force  r/j.  is  replaced  by  where  p  is  the  (h'lisity  in  mass  per  unit  volume 

of  the  shell.  Similar  svvbstitutions  are  made  for  (]$  and  </„.  By  com\)ining  (2.0)  and  (2.7),  one 
arrives  at  the  time-dependent  Donnell-Mushtari  equilibrium  ecpiations  for  a  thin  cylindrical 
sin'll  with  radius  of  curvature  H  and  thickne.ss  h 
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VVe  note  that  the  representation  of  the  <'xternal  loads  as  surface  moments  and  force's  is 
convc'iiient  when  deriving  the  strong  form  of  the  e'epiations  of  nu:)tion.  However,  in  many 
applications  whe're'  it  is  lu'ce'ssary  to  actually  de'te'rmine  e'xpressions  for  these  loads  or  when 
using  the'  we'ak  form  of  the  eepiations.  it  is  advantage'ous  to  re])resent  these'  le)aels  in  te'rms 
of  line  forces  and  moments.  To  accomplish  this,  led  A/j.,  .f/^,  .\(r.  and  .V^  denote  the  external 
re'sultants  acting  on  the  eelge  of  an  infinite'simal  element  which  have'  the  same  orientation  as 
the  internal  re'sultants  elepicte'el  in  Figure  2  (with  units  e)f  me)ment  and  fe)rce  per  unit  length 
of  mieldle  surface).  Fe)rce  and  moment  balancing  can  be'  use'ei  te)  write  the  area  moments  and 
in-i)lane'  fe)rce's  in  te'rms  e)f  tlu'se  line  moments  and  fore'e's,  thus  yiehling 


qr  = 


lUr  =  -■ 


()Nr 
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1  ()i\Io 

R  dO 
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1  <)i\g 

IIW 

()M. 
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We'  point  oeit  that  the  first  exprt'ssion  in  (2.9)  ran  be  obtaiiu'd  from  (2.6)  simply  by 
re'placing  A'j-  by  W-  anel  deleting  iu  the  fir.st  expre'ssieui  eT  (2.6).  .Similar  analysis  leads  te) 
the  e)ther  e'xpre's.sie)ns  in  (2.9).  The'  u.se  of  thf'.se  line  moments  anel  force's  in  (2.8)  is  e'quivale’iit 
te)  inclueling  the  exte'rnal  re'sidtants  directly  wlien  determining  the  e'quations  of  moment  anel 
feerce'  e-epiilibrium  fe)r  an  infinitesimal  shell  eleme'iit  as  ele)ne  in  (2.6)  and  (2.7). 

The'  substitution  of  the  internal  moments  anel  feerce's  in  (2.5)  anel  the  e'xte'rnal  re'sultants 
fre)m  (2.9)  the-n  yielels 
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2R  d.rd0 
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1  ()^m  P  du  1  dv  1  h^  (1  —  p^) 
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Cf  dP  Rd.r  R^  dl)  R^  12  Eh 


Eh  d.r 


(1  -P^)  1  dNg 

Eh  R  dO 


dHh  dMr 


[2.10) 


R^  do^ 


d.r^ 
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where  again,  u,  v  and  w  are  the  axial,  tangential  and  radial  displacements,  respectively  [14]. 


The  constant  Ci  given  bj 


is  the  phase  speed  of  axial  waves  in  the  cylinder  wall.  The  external  line  forces  Nj.  and  Ne 
and  moments  A/j  and  Mq  have  units  of  force  and  momeTit  per  unit  length  of  middle  surface, 
respectively,  and  are  generated  in  our  problem  by  the  activation  of  the  piezoceramic  patches. 
The  load  </„  is  left  as  a  surface  force  since  this  is  the  form  that  it  usually  takes  in  problems 
involving  the  excitation  of  a  shell  through  normal  forces  (an  example  of  a  normal  force  in  this 
form  is  the  pressure  exerted  on  the  shell  due  to  an  exterior  or  interior  acoustic  field). 

We  again  emphasize  that  the  resultant  expressions  in  (2.9)  (and  hence  the  system  (2.10)) 
were  derived  for  an  infinitesimal  element;  hence  certain  modifications  must  be  made  when 
considering  the  global  form  of  the  resultants  and  equations  (as  is  necessary  when  the  resultants 
are  generated  by  a  piezoceramic  patch).  In  certain  cases  (e.g.,  for  certain  types  of  moments 
and  forces),  the  system  (2.10)  agrees  A^ith  the  strong  form  of  the  global  shell  equations.  In 
general,  however,  this  is  not  true,  and  one  must  exercise  extreme  care  in  determining  the  form 
of  the  global  representations  for  the  moments  and  forces. 


Weak  Form  of  the  Donnell-Mushtari  Cylindrical  Shell  Equations 


In  order  to  find  the  weak  form  of  the  shell  equations,  the  kinetic  and  strain  energies  of 
the  shell  are  needed.  By  combining  the  Kirchhoff  shell  hypothesis  with  the  strain  results 
from  classical  elasticity  theory,  it  follows  that  the  strain  energy  stored  in  the  shell  during 
deformation  is  given  by 


1  /■^/2  r2ir  j-e 

::  /  /  +  <70e9  +  <TrOx^)(l  +c//?) 

z  J-h/2  Jo  Jo 


RdxdOdz 


where  the  stains  and  stresses  are  given  in  (2. 1)  and  (2.4),  respectively.  Substitution  and 
integration  (with  zj /?)“'  replaced  by  its  geometric  series  expansion  and  neglecting  powers 
of  c  in  the  integrand  which  are  greater  than  two)  yields 


I  /-^TT  ff  Eh  \ 

"2/0  /o  +  -2(i-//) 


(I  -  ,.^) 


\  2  f  2 

+  T^  (^x  +  K-e)^  ~  2(l  —  //)  (  — T  j  —  £e^e) 


With  the  change  of  variables  .s  =  x/ R,  the  total  strain  energy  can  be  written  as 

1  r2ir  f(/H  fU 

=  /  71 - 7T  +  klsn]  dsdO 

1  Jo  Jo  ( 1  — 
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wlu'n'  h  =  1)^ /{\‘2H^),  IpM  is  tiu*  iiitegraiul  rorn'spoiKliiig  tt)  tlu'  Doiiiicll-M uslit ari  (lirorv 
aiul  Ini-'i  <l('ii()t('s  th('  terms  which  are  retaiiK'ii  to  viehl  th('  Hyriu'.  Fhiggt'  and  Imr've  strain 
eiK'rgy.  Fln'se  two  componenls  arc'  givc'ii  !)y 


On  dv  \  ^  ()u  { ()v  \  1  /  ()v  ()ti 


,  I  /„  )  \  (  ()^n 

+  A.t(V„.)  -■^,1-0  ___  _ 


,  ,  drePw  (Ftr  it  ,  f  ()i 


+  (  1  ~  )  T 


i)ll  i)~U' 


()()  ()s<)0 


1  , ,  ,  / Ou' 


()uifu'  <Pw 

-  2 - 1-  2ir - 1-  ir 

(h  (h^  do^ 


For  siini)ricity  of  presentation,  a  weak  form  of  the  sliell  ecpiations  will  hc'  (h'veioired  using 
the  Donn('ll-M\ishtari  strain  ('xiuession;  a  eorres|)on<rmg  set  of  ecpiations  can  lu'  dc'rived  in  a 
similar  manner  in  the  Byrne.  Fliigge  and  Ltir'ye  case. 

TIh'  kinetic  ('iK'rgy  of  the  shell  is  given  by 


riirouglujiit  this  development,  it  is  assumed  that  the  shell  satislies  slunr  di(ip})raijiii  bound¬ 
ary  conditions  at  =  0,  f;  that  is.  it  is  assumed  that 

V  =  W  =  .V,-  =  Mr  =  U 

at  the  ends.  This  is  done  merely  to  demonstrate  the  c'cpiivalenc  t'  bet  we-en  the  wc'ak  form  which 
follows  and  the  strong  form  already  discussed;  oth.'r  boundary  conditions  can  l)e  treatc'd  with 
similar  argunu'uts.  It  should  be  noted  that  th«'  (onditions  v  =  u’  =  0  at  t  he  ends  are  essc-ntial 
boundary  conditions  and  lu'iice  must  Ix'  enforcc'd  on  the'  chosen  state  s]>a(t'. 

For  an  arbitrary  time'  intc'rval  [/u./i].  considc-r  the  action  intc'gral 

A[i!]=  f\T-V)dl  (2.11) 

Jio 

where  ii  =  [«,?’,  (/;]  is  considc'red  in  the  space'  V'  =  //(Jlf))  x  x  H^(U).  fb'rc'  if  (h'liotc'S 

th('  shell  and  tlu'  subscript  h  denotc's  the  set  of  functions  satisfying  tlu'  c'ssential  boundary 
conditions.  One  then  considers  variations  of  the  form 

it{t.r.0..r)  '  r/,(/)d,(r.fT.r)  ' 

u  =  (7 -f  c<F  =  r{t,r.(),.r)  tf.  .r) 

u’(/.  7'.  0.  .7’)  0.  .r) 

Here  ij  ~  [hi-'/2-hal  '"'d  =  [Oi.d7j,0:i]  are'  chosen  s<7  that 
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i. )  uit,;;-)eV 

ii. )  u{to,  ^  . 


Note  that  tliis  enforces  if  6  [//^(O,  7’)]'\  if(fo)  ~  v{U)  ami  6  €  V  . 

lla?iiilton\s  priticiple  states  that  the  motion  of  tiie  shell  m\ist  give  a  stationary  value  to  the 
action  integral  when  compared  to  variations  in  the  motion,  thus  h’ading  to  t  h<’  r('(|nir('m<'nt 
that  for  all  <?. 

x-' I 


=  0 


f=0 


With  the  definition  (2.11)  for  tlu'  action  integral.  Hamilton’s  principle  leads  to  the  condi  ion 
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.Not('  that  this  must  hold  for  all  arbitrary  intervals  [/o,/i]  and  all  admissible  perturbations. 
Tempt)ral  integration  by  |)arts  in  the  first  integral  in  conjunction  with  tlu'  underlying  condition 


1  1 


that  17(^0)  =  vi^i)  then  yields  the  coi.pled  system  of  equations 
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d^it)  d^4>3  d^w  d'^(j)3  ^  d^w  d^(p:i 
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I  =  0 


The  weak  form  of  the  equations  of  motion  for  the  unforced  shell  is  thus 
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8s89  8s89\  j  J 


for  all  $  E  V.  Again,  the  constant  Cl  =  [^(TTT^]  phase  speed  of  axial  waves  in  the 

cylinder  wall. 

In  terms  of  the  moment  and  force  resultants  (see  (2.5))  and  the  original  axial  variable  .r, 
the  weak  form  is 


r  f  I  = 0 

Jo  Jo  [  8t^  8x  89  J 


--  ^(l>2  .  ^^2] 


r  f  \  +  RN,e~  ^  dxd9  =  0 
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(2.12) 


r2n  rf  {  (Pw  ,  ...  ^’^03  1  (^^3  8^<p2 
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The  derivation  thus  far  has  been  for  the  unforced  shell.  To  include  the  contributions  of 
applied  external  forces  and  moments  which  do  nonconservative  work  on  the  shell,  one  can  ap¬ 
peal  to  an  extended  form  of  Hamilton’s  principle  or  more  formally  include  these  contributions 
directly  in  the  system  (2.12).  Both  techniques  yield  identical  final  equations  and  for  ease  of 
presentation,  we  will  take  the  latter  approach. 

The  inclusion  of  the  applied  line  forces  and  moments  Ng  and  Mj.,  Mg  and  the  surface 
load  (/„  in  the  system  tluMi  yields 


+  +  Ng,^  -  dxdo  =  0 

ox  aO  ox  J 


(  d^w  d^d)-\  I 

/  /  \ 

Jo  2o  ot^  oxJ  R  ou^ 
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-Rq,,4>3  +  RM,^  + 

ox^  R  o6^ 


'dxdO  =  0 


(2.13) 


for  all  0  G  V'  as  the  weak  form  of  the  Donnell-Mushtari  equations  of  motion  for  the  forced 
shell. 

With  the  assumption  of  sufficient  smoothness,  the  weak  solution  in  this  form  is  consistent 
with  the  strong  solution  in  (2.8).  The  vanishing  of  several  of  the  boundary  terms  which  arise 
during  integration  by  parts  is  a  result  of  the  choice  V  =  Rgi^)  x  Rhi^)  ^  Rhi^) 
function  .space  since  the  state  variables  and  test  functions  are  required  to  satisfy  the  essential 
boundary  conditions 

IJ  =z  XU  =  0 

at  O’  =  0,  f. 

We  point  out  that  in  the  weak  form  (2.13),  one  is  not  required  to  differentiate  the  applied 
force  and  moment  resultants  Ng,  and  Mg  as  is  required  in  the  strong  form  (2.10).  This 
proves  to  be  very  beneficial  when  these  terms  are  generated  by  finite  piezoceramic  patches  as 
discussed  in  the  next  section. 


Plate  Equations 

Consider  a  thin  rectangular  plate  whose  edges  lie  along  the  coordinate  lines  x  =  0,  f  and 
^  =  0,  fl.  We  assume  that  the  plate  subjected  to  both  longitudinal  and  transverse  loading 
via  the  surface  forces  and  moments  and  rhj.,xhg.  With  u,x^  and  xu  denoting  the 

displacements  in  the  x,y  and  normal  directions,  respectively,  the  strong  form  of  the  Kirchhoff 
plate  equations  is  given  by 
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where  the  moment  aiul  force  resultants  are 
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The  first  two  equations  in  (2.14)  describe  the  longitudinal  movement  of  the  plate  while  the 
third  equation  describes  the  transverse  motion  of  the  plate. 

To  find  the  weak  form  of  the  equations,  the  vector  u  =  [u.t’.u.’]  containing  the  displace¬ 
ments  in  the  x,y  and  normal  directions  is  considered  in  the  space  V'  =  Hl(Cl)  x  //^(n)  x  //(f(ri) 
where  fl  denotes  the  plate  and  the  std>script  b  denotes  the  set  of  functions  satisfying  essential 
boundary  conditions  for  a  specific  problem.  By  using  analysis  similar  to  that  just  described 
for  cylindric-al  shells,  the  w('ak  foriTi  of  the  equations  of  motion  for  the  plate  can  be  found  to 
be 


0<j) 


0<t>i 
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(2.15) 
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dxdy  =  0 


for  all  (f)  =  [01 ,  (^>2,  d".!]  G  V-  As  in  the  case  of  the  thin  shell,  the  external  line  forces  and 
moments  Nj.,  Ny,  Mj-  and  My  are  related  in  an  infinitesimal  semse  to  the  corresjionding  area 
forces  and  moments  qj-,(iy,ih.y  and  ihr  appearing  in  the  strong  form  of  the  equations  by  the 
relations 
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(2.16) 
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If  the  solution  has  sufficient  smoothness,  integration  by  parts  can  be  used  to  show  that  the 
weak  solution  is  consistent  witli  tlie  strong  solution  in  (2.14). 

Beam  Equations 


The  motion  of  an  undamped  thin  beam  of  length  (  and  width  1  can  be  determined  from 
the  dynamics  of  thin  plate  theory  by  considering  oidy  the  vibrations  in  the  .r-direction  along 
with  the  usual  transverse  vibrations  (in  the  r  direction).  From  (2.14)  this  yields  the  strong 
form  of  the  Euler- Bernoulli  beam  equations 

d^u  dNr 

(2.17) 

d^ui  .  dwy 

^  ~  ^ Ih 


where 


.V, = e4 

ax 


=  -El- 


Note  that  /  =  h^/l2  is  the  moment  of  inertia  for  a  beam  of  width  1. 

A  corresponding  weak  or  variational  form  of  the  equations  can  be  determined  by  choosing 
1/  =  H^(Q)x  U^i^)  for  the  space  of  trial  functions  where  Q  denotes  the  beam  and  the  subscript 
b  again  denotes  the  set  of  functions  which  must  satisfy  the  essential  boundary  conditions. 
Through  either  an  energy  derivation  such  as  that  given  for  the  thin  shell,  or  simply  integration 
by  parts,  one  arrives  at  the  variational  form 


+  A^x^-iV.^|r/.r  =  0  for  all  e//,'(n) 

~  dx  =  0  for  all  ^>3  G 


(2.18) 


of  the  beam  equations.  We  point  out  that  in  this  form,  one  is  not  required  to  differentiate  the 
('xternal  force  or  moment  resultants,  Nj^  and  Mj..  which  proves  to  be  very  useful  when  these 
terms  are  generated  by  the  activation  of  finite  piezoceramic  patches. 


3  Patch  Contributions  to  the  Shell  Equations 


For  a  thin  cylindrical  shell,  the  strong  and  weak  forms  of  the  equations  of  motion  are  given 
by  (2.8)  and  (2.1.4),  respectively.  In  the  case  of  the  weak  form,  it  is  seen  that  the  loads 
can  be  written  in  terms  of  the  line  forces  and  moments  AV-  A(>«  Hiid  Mo  and  the  normal 
surface  load  r/,i,  while  the  strong  form  contains  surface  loads  and  the  derivatives  of  surface 


iiionieiits.  In  the  problem  under  consideration,  these  quantities  result  from  the  activation  of 
piezoceramic  patches  of  thickness  T  which  are  assumed  to  he  perfectly  bonded  to  a  cylindrical 
shell  of  thickness  h  with  midsurface  radius  H  (see  Figure  d).  As  shown  in  Figure  4,  the  patches 
are  assumed  to  he  situated  so  that  their  edges  are  parallel  to  lines  of  constant  x  and  0.  Because 
the  patches  generate  no  shear  strains,  the  exterior  load  q„  is  taken  to  he  q„  —  0.  If  the  weak 
form  (2.13)  is  used,  the  external  line  moments  and  forces  are  simply 

Mx  =  (A/x)p.  <  Me  =  {Mg\,, 

Nr  =  {N,)re  ,  Ne  =  {Ne)j,, 


where  {Nr)pF  and  (Ng)pp  are  the  respective  moments  and  in-plane  forces  which 

are  generated  hy  the  patches.  The  subscript  pc  is  used  to  denote  patch  properties  and  to 
help  differentiate  them  from  shell  properties  which  have  no  subscript.  When  it  is  necessary 
to  differentiate  between  the  two  patches,  the  outer  w'ill  be  denoted  with  a  subscript  pcj  with 
a  subscript  pe-z  being  used  to  denote  the  inner  patch. 

However,  if  one  is  using  the  strong  form  (2.8)  of  the  equations  of  motion  with  piezoceramic 
actuators,  the  surface  moments  and  forces  to  be  used  in  (2.8)  are  given  by 


irir  =  —  ■ 


1  d(Mg)pF 

R  30 


mg  = 


d{M.) 


pe 


dx 


diNr)pe 


(]e  - 


d[Ne) 


(3.2) 


pf 


For  a  patch  with  uniform  thickness  and  bounding  values  X),  and  Oi.  the  presence  of  the 
indicator  function 

■  1  ,  j-  <  (.ri  +  x-i)!'! 

0  ,  .r  =  (.r, +.r2)/2  (3.3) 

-1  ,  .r  >  (ji -f  j-2)/2 


derives  from  the  fact  that  the  forces  generated  by  the  patch  in  the  j-direction  are  antisym¬ 
metric  (equal  in  magnitude  but  opposite  in  signj  about  the  line  x  =  (ji  -|-  x-i)!'!.  The  same 
holds  true  for  the  forces  in  the  ^-direction  with  S\,2[0)  being  defined  in  an  analogous  manner 
to  .S',  2(.r)  in  (3.3). 

We  point  out  that  the  differences  between  the  external  surface  force  expressions  in  (2.9) 
and  (3.2)  are  due  to  the  fact  that  the  former  were  derived  for  an  infinitesimal  element  whereas 
the  latter  are  global  expressions  which  preserve  the  overall  signs  of  the  forces  generated  by  the 
patches  as  well  as  reflect  the  discontinuities  d»ie  to  changes  in  sign.  These  differences  result 
from  the  property  that  the  sense  of  the  forces  is  highly  dependent  on  the  specified  location 
of  the  axis  origin  on  the  neutral  surface.  Hence  the  direction  of  forces  throughout  the  patch 
fliffers  in  .some  locations  from  tho.se  okserved  in  the  infinitesimal  element  thus  necessitating 
the  inclusion  of  the  indicator  functions  in  (3.2). 

Unlike  the  forces,  the  action  of  the  moments  is  specified  with  respect  to  a  fixed  point  on 
the  tK'utral  surface  (the  point  0  for  the  element  in  Figure  2  or  a  point  on  the  left  edge  of  the 
shell  in  hdgure  1).  As  long  as  the  orientation  of  the  infinitesimal  element  and  full  shell  with 
patches  are  the  same,  the  line  moments  derived  for  the  infinitesimal  element  will  be  consistent 
with  those  of  the  full  structure.  Thus  the  expressions  for  the  general  infinitesimal  moments 
in  (2.9)  nee<l  no  modifications  when  clescribing  the  surface  moments  generated  by  the  patches 
as  given  in  (3.2). 


If) 


Figure  4.  Piezoceramic  Patch  Placement. 


In  order  to  determine  {Mx)pe,{^9)peA^x)pe  and  {Ng)pe  and  hence  the  loads  on  the  shell, 
it  is  useful  to  write  them  in  terms  of  the  normal  strains  ex,£e  and  midsurface  changes  in 
curvature  Kx,Kg  of  the  middle  surface  {z  =  0)  of  the  cylindrical  shell.  That  is,  we  want  to 
express  the  patch  moments  and  forces  in  terms  of  the  reference  surface  characteristics  of  the 
cylindrical  shell. 

We  emphasize  that  due  to  the  presence  and  activating  nature  of  the  patches  when  a  voltage 
is  applied,  the  normal  strains  and  changes  in  curvature  are  no  longer  given  by  the  expressions  in 
(2.3)  which  were  derived  for  a  homogeneous  thin  cylindrical  shell.  At  this  point,  £x,£e-,  and 
K$  are  considered  to  be  unknown  and  are  determined  by  formulating  stress-strain  relations  in 
the  patch  and  shell  followed  by  the  balancing  of  moments  and  forces  in  the  combined  structure. 
In  this  manner,  expressions  for  these  midsurface  characteristics  (and  hence  the  resulting  patch 
moments  and  forces)  can  be  found  in  terms  of  the  material  properties  of  the  shell  and  patch, 
the  radius  of  curvature  of  the  shell,  and  the  voltages  being  applied  to  the  patches. 


Stress-Strain  Relations  in  the  Patch 


Prom  (2.1),  the  infinitesinially  exact  strain  relationships  for  a  cyliiKlrical  shell  vvitli  mid- 
snrface  radius  R  are  given  by 

fx  =  (tr  +  =Kr) 

-  TTTJii 

where  Cj.  and  eg  are  the  normal  strains  at  an  arbitrary  point  on  the  cylindrical  shell.  If 
the  patches  and  cylindrical  shell  are  thin  in  comparison  with  the  radius  of  curvature  of  the 
shell,  then  it  is  reasonable  to  assume  that  the  relationship  (2.4)  is  maintained  throughout  tlu' 
combined  thickness  h  +  2T  as  shown  in  Figure  3  {.see  also  [12]).  Hence  we  will  tak(’ 

(^v)pf  —  ('ix  d"  ~^x) 

{('e)pe  —  j  {£e  +  =Kfl) 

where  (fj-)pp  and  (co);,^  are  the  normal  strains  at  aji  arbitrary  |)oint  on  tin'  patch.  Note  that 
this  assumption  implies  that  the  strains  at  the  interface  an'  continuous  and  that  the  (('iitc'is 
for  the  radii  of  curvature  for  the  shell  and  patch  are  concurrent.  As  seen  in  (2.4)  and  (.2.5). 
the  tangential  strain  distribution  in  the  shell  and  patch  is  in  general  nonlinear  in  r  (see  also 
Figure  2). 

VVe  point  out  that  the  model  at  this  point  differs  from  the  flat  ])late  piezoceramic  coupling 
model  [15]  both  in  the  presence  of  the  term  z/R  in  the  tangential  straiti  expn'ssion  of  (2.5) 
and  in  the  fact  that  the  model  retains  the  coupling  between  the  normal  midsurface  strains  and 
the  changes  in  curvature  (this  is  analogous  to  simultaneously  considering  both  longitudinal 
and  transverse  vibrations  in  a  plate).  Although  the  ratio  z/R  is  neglected  when  deriving 
the  Donnell-Mushtari  model  (see  (2.2)),  we  retain  it  here  so  that  curvature  elfects  are  fully 
included  in  the  coupling  between  the  patch  and  shell.  The  retention  of  this  term  also  ensures 
that  the  patch  ititeraction  model  can  be  directly  applied  to  higher  order  shell  models  without 
necessitating  changes  to  accommodate  the  greater  accuracy. 

From  the  constitutive  relations  in  (2.4),  it  can  be  seen  that  the  stress  distributions  within 
the  cylindrical  shell  are 

E 

^x  =  - T  (fx  +  t'ce) 

‘  -  "  Cl.C) 

=  - - 7  {eg  +  //e^)  . 

1  — 

The  stress  distribution  in  the  patches  will  contain  contributions  from  both  the  free  piezo- 
ceramic  actiiator  strain  and  the  strain  distribution  in  (2.5).  At  this  point  we  assume  that 
when  voltage  is  applied  and  the  patch  is  activated,  in  accordance  with  basic  shell  theory, 
equal  strains  are  induced  in  the  x  and  0  directions  and  the  radius  of  curvature  is  not  changed 
in  either  direction.  Patches  satisfying  this  assumption  could  be  made,  for  example,  by  taking 
a  portion  of  a  thin-walled  tubular  piezoceramic  element.  For  the  outer  patch,  the  Tiiagnitude 
of  th('  induced  free  strains  is  then  taken  to  be 

(  —  (f  x)pri  =  (f  @)pr|  =  ^  * 


IS 


where  (i^  is  a  piezoteraiiiie  strain  constant  ami  I'l  is  the  ai)plie(l  voltap,<'.  We  point  out 
that  when  a  voltage  is  applied  to  a  patch  with  edge  coordinates  .r\,.r2.0]  and  O2.  tin'  point 
(.vj))  =  ((.rj  +  d-2)/‘2,  H{0]  02)12)  will  not  move  wln'reas  tlu'  axially  symmetric  i)oints  on 

('ither  side  will  move  an  equal  amount  in  opi)osite  din'ctions.  I’his  motivatc's  the  use  of  the 
indicator  functions  at  various  points  throughout  the  d<'velo|)m('nt 

Assuming  that  the  two  i)atches  have  the  same  Young's  modulus,  /v^,, .  and  Poisson  s  ratio, 
//p, .  the  stress  distrilnition  in  the  outer  patch  is  given  by 

“j  ^  j-  d”  9  (Id"  ^''pi  )^  pf\  ) 

K,  ' 

{^9) pf  i  "j  ^  ^  d“  t'^pt  ^  j-  (1  d"  ;'f  1  ) 

with  tlu'  m\gativ('  signs  resulting  from  conscuvat ion  of  forces.  Similar  c'xpressions  arc'  usc'd  for 
tin'  induced  free  strain  and  stress  distribution  in  the  innc'r  ])atch.  By  c-omparing  (d.b)  and 
(d.T),  it  can  be  noted  that  a  stress  discontinuity  occurs  at  tlu'  intc'ifacc'  dm'  to  tlu'  dilfc'rc'iit 
material  properties  of  th<'  shell  and  patch. 


Moments  and  Forces  in  the  Patches 


By  integrating  the  stressc's  over  the  face*  of  a  fundamental  c'leim'iit.  it  follows  that  the' 
moment  and  force  residtants  for  the  patches  can  be  expre.ssed  as 


fii/2+r 

=  /  ( 

Jk/2 

rh/‘2+T 

{^^o)pf\  —  /  )^$)p(\~0z 

J  It./ 2 

flt/2-pT  /  ~  \ 

(c'VW.  Mpr,  +  = 

fii/2vr 

Jk/2 


r  I  pf  I 
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0>f  I  dc 


{‘^h)pf2  -  (’  ■*■  /y) 

r-hl2 

(■^Jo)p<2  ~  I  (^f))pf2~'^~ 

J  —  hl'2  —  T 

r-hl2  / 

(  1  d- 

J-h/2-T  \ 

f-k/2 

(Ao)j„2  /  (c7())pf;r/_ 

J~h/2-r 


(d.s) 


with  units  of  moment  per  unit  length  and  force  p('r  unit  length,  respc'ctivc'ly.  TIk'  c'xplicit 
depc'iulence  of  the  patch’s  moment  and  force  resultants  on  thi'  sin'll  miclsurface  characterist  ic  s 
and  K0  can  be  seen  by  combining  (3.r)),  (3.7)  and  (3.(S)  (we  again  point  out  that  tlu' 
midsurface  strains  £j.,ee  and  curvature  changes  Kj-.ko  are  unknown  ami  will  In*  detc'rmined  by 
balancing  moments  and  forces  in  the  combined  structiire). 


Patch  Loadings 


It  should  b('  noted  that  throughout  this  d<'v<>lopment .  t'dgc'  t'lfects  dm'  to  the'  pate  hc's  have' 
been  ignored  and  thus  the  exi)r('ssions  in  (3.S)  apply  to  patchc's  covc'iing  the'  full  c  ircumfc'n'ncc' 
of  the  shell  and  having  infinite  axial  length.  The  etpiations  can  be'  moclilie'd  for  linite'  |)atc  lie's 


in  the  following  manner.  For  a  patch  with  hounding  values  .rj,  X2,  0i  and  O2  as  shown  in 
Figtire  4,  the  total  line  moments  and  forces  are 

(M,),,,  =  +  (Me)peJ  [H,(X)  -  //-.(x)]  [H,{9)  -  //,(0)] 

{Me)re  =  [(A/.)pe,  +  (A/.VJ[//,(.r)  -  W,(.r)]  (//,(<?)  -  //-.(f?)] 

(;f9) 

(W.)pe  =  [{KU  +  [//.(■^)  -  -  H2{0)]S^a(x)S,,2(O) 

(NeU  =  +  (^e)peJ[//,(^)-  //2(j-)][//i(0)  -  H2mS,a(x)Su2{0) 

where  H  is  the  Heaviside  function  and  Hi{x)  =  H{x  —  Xi)  ,  f  =  1,2,  with  a  similar  definition  in 
9.  The  indicator  functions  .S'i,2(^)  ^md  (see  (3.3))  again  derive  from  the  property  that 

for  homogeneous  patches  having  uniform  thickness,  opposite  but  equal  strains  are  generated 
about  the  point  {x,9)  =  ((x|  -f  J^2)/“>  /?(^i  +  ^2)/^)  hi  the  two  coordinate  directions. 

The  combination  of  the  expressions  (3.8)  and  (3.9)  yields  the  patch  moments  and  forces 
('Afx)pe,  ( A/0)pf,  ( A^i.)pe  and  (A^fl)pf  in  terms  of  the  middle  surface  characteristics  of  the  the 
cylindrical  shell.  Integrating  the  expressions  in  (3.8)  is  somewhat  cumbersome  however,  and 
the  procedure  can  be  facilitated  by  determining  the  patch  moments  and  forces  in  terms  of  the 
resultants  of  the  forced  shell.  To  accomplish  this,  force  and  moment  balancing  is  employed. 


Determination  of  the  Patch  Moments  and  Forces 


The  application  of  moment  equilibrium  about  the  center  of  the  shell  yields  the  two  condi¬ 
tions 

A/x  +  (Mx)pe,  +(A/x)  pe2  ~ 

Me  ^  {MeU  +  (Mg),,,  =  0 

where  Mx  and  Mg  are  shell  moments.  Similarly,  force  equilibrium  in  the  x  and  9  directions 
yields 

A/x  +  (A^i)pe,  +  (A^r)pe2  ~  ^ 

+  i^e)pei  +  {^e)pe2  =  0  ■ 


Thus  the  total  patch  resultants  can  be  expressed  as 


(A/x)pr|  T  {Mx)pe2  —  Mx  —  j  li/-  0 

fhl2 

{Mg)pej  +  {Mg)px2  =  -Mg  =  -  /  agzdz 

J~hl2 

(^T)pn  +  (A^x)pf2  =  =  ~  j 

rhl2 

{^0)pe,  +  i^0)pc2  =  -Ng  =  -  /  agdz 

J-h/2 


(3.11) 
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The  two  tangential  expressions  are  approximate  in  tl)e  sense  that  the  terms  (1  +  z/ R)~^  are 
replaced  by  the  truncated  geometric  series  X—zjR  before  integration  (this  is  the  same  strategy 
which  is  used  when  determining  tlie  moment  and  force  resultants  in  the  Byrne-Fliigge-Lur'ye 
general  shell  theory).  The  patch  resultants  in  (3.12)  are  then  used  in  (3.9)  to  determine  the 
total  line  forces  and  moments  generated  by  the  finite  patches.  Because  these  resultants  are 
functions  of  the  material  properties  as  well  as  the  midsurface  characteristics,  they  can  be 
easily  constructed  once  tx,ee,Kr  and  have  been  determined.  This  is  again  accomplished 
by  moment  and  force  balancing. 


Determination  of  the  Midsurface  Characteristics 


In  terms  of  the  stresses,  the  moment  and  force  equilibrium  equations  (3. 1 0)  and  (3. 1 1 )  can 
be  written  as 


The  integrals  appearing  in  (3.13)  and  (3.14)  are  explicitly  evaluated  in  the  [1]. 
After  collecting  terms,  this  yields  the  4x4  linear  system 


[As/iW/  +  A]  +  A-i]  e  —  epc/i  +  f  -pe2f'i  (3.15) 

•  T  rr’ 

in  the  unknowns  c  =  (Sx,  £0,  Kg)  .  The  shell  contributions  in  Asheit,  the  outer  patch 
contributions  in  A\  and  /i,  and  the  inner  patch  terms  in  A-i  and  are 
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i  =  1,  •  •  • ,-'},  are  of  order  at  most  three  in  h. 


Note  that  the  roeffirients  a,  , 


We  peiiil  out  that,  /li. /Vj, /i  and  J2  dejx'iid  on  inatt'rial  inoperties  of  tlie  patclu's 

and  sliell  (tliickiu'ss.  elaslir  j^ropert irs,  and  [’oisson  ratios),  tlie  radius  of  eurvatnrc'  of  the 
slu'll,  and  the  voltage  Ix'ing  aiiphu'd  to  tlu'  patelu's  (recall  that  =  (I:n\\j"l'  where'  V)  is  the 
a|)i)li('d  voltage'  into  the'  e)nt('r  patch  with  a  similar  ele'finitie)n  for  e,,,., ).  d'he  above'  fejnnulat ie)n 
iseelate's  the'  cont rihut ions  ehu'  te)  the'  inelivielual  patclu's  an<l  is  nse'fnl  if  eeiu'  wants  te>  activate' 


Algorithm  for  Determining  the  Shell/Patch  Interactions 

The'  ste'ps  whieh  are  nee'e'ssary  for  solving  for  the  shell  loads  dne'  te)  the'  activation  of  the' 
patclu's  can  he  snmniarizeel  as  follows. 

(1)  Se't  np  the  1  x  1  syste'm  Ac  =  /  and  solve  for  c  =  Hj..  koY  which  contains  the 

midsurfae'e'  strains  and  change's  in  curvature'. 

(2)  De'te'rniine'  the  line  moments  anel  forces  which  are'  generateel  by  the'  inelivielual  patclu's  as 

set  np  in  (d.  12). 

(3)  The  ce)rresi)e)ndiiig  e  e)mbined  resultants  for  finite'  jiatche's  are'  give'ii  by  (3.!)). 

(4)  riie'  re'snitants  from  (3.9)  are  dire'ctly  substit  ute'el  into  (2.13)  as  t  he  le)ael  on  the'  she'll  if  the' 

we'ak  she'll  eepiations  are  being  use'd  (re'call  that  in  this  case,  e/„  =  0  anel  .V,.  —  (.Vj.)^„, 
.\f)  =  (Ao'jjir,  Mx  =  ('V/,);.f, .su  HI  Ilia  lizeul  in  (3.1)).  For  the  stre)ng  feuni  e)f 
t  he'  e'finatieens  of  motion,  the'  elerivative  e'X|)re'ssie)ns  in  (3.2)  are  fornu'el  anel  snbst  it  ute'el 
inte)  (2.8)  as  the'  external  lejael. 

We'  pe)int  eeut  that  the'  substitution  e)f  the'  pateh  meune'iits  anel  fence's  inte)  the'  stre)ng  feerm 
of  the'  she'll  eepiatie)ns  re'sidts  in  ecne’  elerivative' e>f  the'  lle'avisiele  anel  ineliea1e)r  fune  tieens  fe)r 
the'  fe)ree's  anel  twe)  eh’rivatives  e)f  the'  Meavisieh'  fune'tie)n  for  the  monu'nts  wlu'ie'as  lu)  siu  h 
eliffe'ientiat ie)n  is  re'epiire'el  in  the'  we'ak  form  (the  elerivatives  are'  t ransferre'el  eenle)  the'  te'st 
functions  anel  eeiu'  simply  inte'grate’s  over  the*  re'gie)n  ce)ve're'd  by  the  leatelu's).  Ihis  is  eeiu' 
motivatie)!!  fe)r  using  the'  we'ak  fe)rm  of  the'  she'll  e'e|uations  in  many  applie  at  ieens. 


4  Patch  Contributions  to  Plate  and  Beam  Equations 

Analysis  similar  to  that  used  for  the  thin  cylindrical  sheik  can  be  used  to  determine  the  forces 
and  moments  which  are  due  to  the  activation  of  piezoceramic  patches  which  have  been  bonded 
to  a  flat  plate  or  beam. 

Plate/Patch  Interactions 

The  patch  interactions  with  a  flat  plate  can  be  determined  in  a  manner  similar  to  that 
used  in  the  study  of  the  interactions  between  a  thin  cylindrical  shell  and  a  pair  of  piezoceramic 
patches  as  discussed  in  the  last  section.  Direct  force  and  moment  balancing  leads  to  the  4x4 
system  listed  under  Method  1,  and  this  system  can  then  be  solved  for  the  unknowns 
and  Ky.  By  then  substituting  these  values  into  resultant  expressions  similar  to  those  in  (3.12), 
one  obtains  the  forces  and  moments  generated  by  the  patch.  This  procedure  can  be  simplified 
however,  by  noting  that  the  strains  in  the  x  and  y  directions  of  a  homogeneous  flat  plate  are 
ecpial  when  equal  free  strains  are  generated  by  the  patch  (see  [6]).  Hence  tw'o  of  the  variables 
can  be  eliminated  which  leads  to  the  more  easily  solved  2  x  2  system  given  under  Method  2. 
It  should  be  noted  that  the  two  methods  yield  the  same  final  force  and  moment  resultants. 

Method  / 

Force  and  moment  balancing  similar  to  that  used  in  the  study  of  the  patch /shell  interac¬ 
tions  yields  the  system 

[-Opiate  +  +  A-2\  e  =  Cpe, /i  +  Cpej/i 

'T 

where  r  =  (Sj-.Sy,  Kj.,  Ky)  .  The  matrices  and  vectors  containing  contributions  due  to  the  plate 
and  two  patches  are 
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with 


Note  that  these  iiiatriees  and  vectors  are  identical  to  those  of  (3.15)  if  one  takes  /?  — >  oc  in 
the  latter  e.xpressions.  Once  £j-,£y,Kx  and  Ky  have  been  determined,  the  moment  and  force 
resnltants  can  be  found  in  a  manner  similar  to  that  used  for  shells. 


Method  2 


Here  we  take  advantage  of  the  fact  that  for  a  homogeneous  plate,  the  strains  in  the  .r  and  y 
directions  will  be  equal  when  generated  by  equal  free  (unconstrained)  strains  from  the  patch. 
Hence  we  take 

f  =  fr  =  =  £  +  K  = 

as  the  strain  distribution  in  both  the  plate  and  the  patch.  This  yiehis  the  stn'sses 


<T  =  (Tx  =  o-y  — 


E 

1  -  // 


in  the  plate  and 

E, 

{^)pe\  {^T)pei  —  {^y)pri  —  (f’  ^ pi"!  ) 

1  l■'pr 

E, 

i'^)p<’2  ~  (^r)pf2  ~  ~  "i  ~  0'''2  ) 

in  the  patch.  The  force  and  moment  resultants  for  the  patches  can  be  found  either  by  int('- 
grating  the  stresses  over  the  patches  or  by  using  force  and  moment  balancing  to  express  (hem 
in  terms  of  the  resultants  for  the  forced  plate.  As  was  done  in  the  shell  analysis,  we  will  take 
the  latter  approach  since  it  yields  simpler  expressions  for  the  external  resultants.  Force  and 
moment  balancing  in  conjunction  with  integration  of  the  forced  plate  resultants  then  yields 


-  EIE 

(tV/j-)p(.,  T  {Mx)pF2  ~  (^^y)pc\  T  (^A/)pf2  -  J  I 

(tVx)pf,  +  {•^x)pc2  ~  (^!/)pr,  +  {^y)pf’2  ~  j  _  • 


(b!) 


As  expected,  these  relations  agree  with  those  in  (3.12)  for  the  forced  shell  with  the  exception 
of  the  0(\/R)  terms  in  the  latter  ca.se  which  are  due  to  the  curvature. 

The  total  resultants  generated  by  a  pair  of  patches  with  edges  parallel  to  lines  of  constant 
.r  and  y  can  be  determined  in  a  manner  similar  to  that  used  with  the  shells.  For  a  patch  with 
bounding  values  .ri,.r2,i/i  and  .1/2,  the  total  resnltants  are 


(.\/x)p,  =  (A4)p,  =  [(A/x)pn  +{M.r)pr2]  [//.(.r)  -  ll2(-r)]  [//,(?/)  -  H2(y)] 

(  .Vr)p.  =  {^y)pr  =  [{Nx)p,  ,  +  ( AT)p.,]  [//.(.r)  -  //2(.r)]  [//,(,/)  -  UAll)]  .S,2(.r)>’,.d.v) 


(1.2) 
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vvlu'H'  again,  H,{x)  =  /7(.r  —  .r,)  ,  i  =  1,2,  .‘'i.2(.r)  (U'nott's  the  iiulicator  funrtiun  (iescriht'd  in 
(.'{..'1),  and  //,(//)  and  .‘>'1,2(2/)  <l<dined  in  an  analogous  manner. 

As  before,  t  and  k  must  he  found  in  order  to  dc'termiiK’  the  resultants  in  (-1.2).  d’his  is 
aecomplished  via  moment  and  force  balancing  which  then  yi('lds  the  system 

[/'Iji/n/f  +  A|  +  /I2]  ('  —  f  pnf\  +  (  pfifi  (  b-^) 

r  r 

where'  (  =  ■  The  component  matrices  and  vectors  art' 
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rile  subscripts  1  and  2  again  refer  to  the  outer  and  inner  patch  contributions,  respectively, 
and  the  constants  Ei,  /?2»'^2  and  ^i.3  art*  given  by 

^  fie(l  -ly)  ^  Eh(l  -  ,y) 

'■  l2E,„(l-r)  •  E,„(l-i.) 


=  s  U-tr  -/c  =  jer  +  y-r  + -r 


«:>=(  '(5  +  ^)  =5r(/,  +  7’). 


Algorithm  for  Determining  the  Plate/Patch  Interactions 

As  in  the  case  of  the  shell,  the  steps  necessary  for  determining  the  plate  loads  which  are 
tint'  to  th('  patcht's  can  bt'  summarized  in  a  simph' algorithm. 

(1)  Set  up  and  solve  tlu'  2 x2  system  A(  =  /  in  (4.3)  where  r  =  (i,  contains  the  midsurfaci' 

strain  and  change  in  curvature. 

(2)  Determine'  the'  combined  re'snltauts  for  the  finite  patches  through  (1.2)  in  conjunction 

with  (-1.1 ). 


(3)  Once  (h'termined,  tlie  resultants  from  (4.2)  can  be  substituted  directly  into  the  weak  form 
of  the  plate  equations  (2.15)  as  the  load  on  the  system  (with  f/„  =  0  and  A^  =  (A’j.),,f. 

Mj-  =  My  =  (My)pe)-  If  the  strong  form  of  the  plate  equations  is 

Ix'ing  us('d,  th('  surface'  heads  can  be'  eh'te'rmine'd  via  the'  expre'ssieens 


-■S.2(.r).S,2(.f/) 


7llj. 


A{My)p, 

(hi 


<h  =  -■'"’i.2(.r)'S,2(.v) 


An 


anel  the'se'  latte'r  value's  eaii  be'  substitute'el  intee  the'  e'rpiilibrium  eepiatieens  (2.14). 


As  in  the  case  of  the  shells,  the  use  of  the  strong  form  results  in  up  to  two  derivatives  of  the 
Heaviside  function  whereas  the  use  of  the  weak  form  alleviates  this  problem  by  transferring 
the  derivatives  onto  the  test  functions. 

It  should  be  noted  that  the  v^oltage  choice  Cpe  =  Cj,,,  =  causes  pure  extension 

(patch  pairs  excited  “in  phase”)  in  the  plate  while  pure  bending  occurs  with  the  choice' 
Cpp  =  — Cpp,  =  Cppj  (“out  of  phase”  excitation). 


Special  Case:  Cpe  =  — Cpr,  =  fpe2 

In  this  case,  the  constants  £  and  k  have  the  values 


t  =  0  , 

which  leads  to  the  total  patch  i7ioments 
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This  line  moment  expression  is  equivalent  to  the  relation 

Pz{2  +  p,) 


A/x  =  My  = 


1  +  d/);(3  +  pi  +  :]p,) 


l,^^cp,[lh{.r)-  H2{j')][HAy)-  Ihiy)] 


with  p;  —  T/h,  h  =  h/'2  which  was  obtained  by  Kim  and  .Jones  [12]  in  their  developnu'iit  of  a 
model  for  the  bending  interactions  between  a  flat  plate  and  a  piezoelect ic  actuator.  In  their 
work,  they  consider  a  patch  configuration  which  excites  pure  bending  in  the  plate  and  then 
determine  the  effective  patch  moment  by  first  isolating  the  interface  stress  of  the  systi-m. 


Beam/Patch  Interactions 


The  patch  contributions  to  the  dynamics  of  a  thin  beam  can  be  determined  diri'ctly  from 
the  plate/patch  interaction  model  if  one  considers  only  vibrations  in  the  .r-direction  along 
with  the  usual  transverse  vibrations.  The  .system  for  the  beam/patch  configuration  is  then 
identical  to  that  found  in  (4.3)  with  the  constants  E\  and  E2  now  given  by 


E,  = 


E 


VIE, 


pc 


E2 


Once  e  and  k  have  been  determined,  the  force  and  moment  resultants  for  tlx'  patch  an- 
expressed  in  terms  of  those  of  tlu*  forced  beam  and  are  given  by 


(A/x)pf,  +  (Mr)pF2  —  —EIk 

(A^x)p<-,  +  ('Vx)pr2  =  —Ehe 
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where  /  =  h'^  j  VI  is  t  lu'  nioiiieut  of  inert  ia  fur  the  beam.  For  patches  with  bounding  values  ,ri 
and  x-i-  the  effective  moments  and  forces  are 


(A/,.)p,  =  [(A/,.)„,  +  (A/,.)pc,]  [//(.r  -  .r,)  -  H[x  -  ,c,)] 

[^x)pr  =  [{l^T)pe,  +  (Avlprj]  [//(•'•  “  J"l)  “  H  {x  -  .r2  )]  -S/il-r) 

which  can  then  be  substituted  directly  into  the  weak  equations  (2.18)  as  loads  on  the  beam 
(with  r/„  =  0  and  AV  =  (AV)^(,  A/j.  =  {Mx)pc)-  In  order  to  determine  the  patch  loads  for  the 
strong  form  of  the  beam  equations,  tlie  corresimnding  surface  moments  and  forces  are  found 
via  the  relationships 


i^h)pe 


^(A(r),,f 

dx 


d{Mx)p, 

dx 


and  these  latter  values  are  used  in  (2.17).  We  again  point  out  that  this  results  in  the  need  to 
differentiate  the  Heaviside  function  (once  for  the  force  and  twice  for  the  moment)  whereas  this 
problem  is  avoided  in  the  weak  formulation  since  the  derivatives  are  transferred  onto  the  test 
functions.  In  fact,  the  effect  o:  the  Heaviside  functions  in  the  latter  case  is  to  simply  restrict 
the  integrals  to  the  region  covered  by  the  patches. 


Special  Case:  Top  Patch  Activation 

Consider  the  problem  of  a  b('am  having  only  a  lop  activating  ])a1ch. 
case  is 
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which  implies  tliat  the  midsurface  constants  are  then  given  by 
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These  expressions  for  e  and  k  are  the  same  as  tho.se  found  by  Cibbs  and  Fuller  [10]  when  they 
wt ’•e  investigating  the  moments  and  forces  generated  by  a  single  patch  which  was  bonded  to 
a  thin  beam  (these  expressions  are  equivalent  to  their  (7)  and  (8)  once  the  latter  have  bet'ii 
simplified  and  h/2  has  been  substituted  for  h).  The  moment  and  force  resultants  (AC),,,,  and 
(Afj.)p(|  can  then  be  found  by  substituting  £  and  k  into  (  l.  I)  (with  {Mr)p(^  =  (A’j.),>f2  =  0  since 
there  is  no  bottom  [latch).  We  note  that  the  resultants  in  this  case  are  ecpiivalent  to  those  in 
[10]  although  the  forms  and  signs  differ  slightly  due  to  a  slight  difference  in  the  formulation 
of  the  underlying  beam  equations. 


5  Conclusion 


I n  t  his  work,  gc’ticral  tiiuch'Is  (h'scrihiiig  t  ho  iiit<'rarf  ioii.s  hoi  w('('ii  a  i>air  of  pi('Zo<'oraiiiio  pat  chos 
and  ail  uiidotiying  (dast  io  stnioliua'  havi'  Ix'oii  |)r<'sonf  od.  W’ldlo  iho  prosoiit  at  ion  is  for  olastii 
siihst iiiot uri's  ooiisisting  of  a  thin  oyliinlrioat  shot!.  plal<'  and  Ixvini.  tho  toohniipn's  discnssod 
tor  dot  orniining  t  lu'  niotnont  s  and  hircc's  goiu'iatod  liy  tho  pat  chos  can  ho  diriat  ly  ('xtondod  to 
mori'  comph'X  strnrtnros  and  gi'ornot  ri('s. 

In  lli('  ('as('  of  th('  sin'll,  tho  |)atohos  ar*'  assunu'd  to  !«'  onrvod  and  tin'  ooniiling  hot  woon 
tho  in-plano  strains  and  the  In'iiding.  whioli  is  dm'  to  tin'  nirvatun'.  is  ri'taiin'd.  By  using  h^rci' 
and  inonn'iit  halanring  to  dc'torniino  tho  midsurfao('  strains  and  changos  in  ourvatuix'  of  tin' 
coniiiiin'd  strnolun',  ('xprossions  for  (In'  [)atrh  inomont  and  Ajiti'  rosnltants  can  ho  d('\ <'l(;pod. 
In  this  nianin'r,  tin'  loading  dno  to  an  aotuating  ])air  of  ])atrhos  can  In'  ('xprc'ssc'd  in  tt'inis 
of  tin'  niatorial  proj)orti('s  of  tin'  slioll  and  patolios  (thirknoss.  olastio  propr'rtii's  and  I’oisstni 
ratios),  tin'  radius  of  ourvaturo  of  tho  slu'll,  and  tho  voltago  hoing  appliod  to  ('arh  of  tho 
patchos.  This  pro\’irlo,s  a  sin'il/pafoh  intoraotion  modo!  which  rotains  tlio  ciirvatiirr  ('IFcrts 
as  \v('ll  as  admits  poti'ntially  dilForont  voltagos  into  (ho  two  ])atrhos.  Wo  point  out  that  tin' 
gt'in'ial  torhnir|uos  ns<'d  for  dotormining  this  cylindrical  sholl/patch  modol  can  also  ho  iisi'd 
to  d('tormin('  tho  intoractions  hi'twoon  pairs  of  piozororamic  patches  ainl  nioro  goin'ral  shells 
(for  ('xam])l('.  in  tin'  case  of  a  spin'rical  shell,  one  woidd  retain  the  cnroatnn'  ofrocts  in  both 
coordinat('  directions). 

Tho  tochni(|uos  for  dc'tormining  tho  patch  intoractions  with  a  cylindrical  shell  wc'ix'  then 
ns('d  to  devc'lop  gc'in’ral  inti'iaction  models  for  |)alches  which  arc'  hoinh'd  to  thin  Hat  jilatc's 
and  In'ams.  .As  in  tho  sin'll  case',  tin'  models  are  snfficic'iit ly  general  to  allow  for  |)ot('ntial!y 
dilfering  patch  voltage's  which  implies  that  they  can  he  nsc'd  for  controlling  systc'm  dynamics 
win'll  both  th'Xiiral  and  extc'nsional  vibrations  arc'  prc'sc'iit.  do  compare  with  existing  aiialysc's. 
tin'  |)late/i)atch  interaction  modol  is  shown  to  ho  c'cptivalc'tit  to  that  of  [12]  in  the'  s|)c'cial  case' 
whc'ii  pure  hc'iiding  motion  is  c'xcitc'd.  .Also,  the'  hc'am/patch  moclc'l  is  c'cpiivalc'iit  to  that  of 
[10]  when  thc'rc'  is  one'  actuating  patch,  llc'ncc'  the*  hc'am  and  jclatc'  intc'raction  modc'ls  arc' 
c'onsistc'iit  with  existing  thc'ories  in  the'  spc'cial  case's  prc'vionslv  oxaminc'd  while  also  allowing 
for  mcnc'  gc'iic'ral  struc  t  ure'/ patc  h  intc'ractions  which  can  arise'  in  more'  complc'x  ap|)lic'a1  ions 
(for  c'xamplo.  coui)lc'cl  systems). 

For  c'ach  of  the'  she'll,  plate  and  hc'am  intc'raction  modc'ls.  the'  contributions  of  tho  patchc's 
are  carefully  doscrihed  in  both  the*  strong  and  wc'ak  forms  of  the'  t imc'-dc'pc'iidont  structural 
c'C|uations  of  motion.  This  jmovidc's  models  which  can  Ice-  iisc'd  in  a  varic'ty  of  apiilicat ions 
including  nimu'iical  simulations,  paramc'tor  iclc'iit ilicat ion.  and  control  schc'mc's.  In  c'ach  of 
thc'sc'  aiiplicat ions,  tho  models  arc'  siifliciontly  gc'iic'ial  to  provide'  for  a  varic'ty  of  approximation 
tc'chnic|uc's  including  modal.  si)c’ctral,  spline'  and  linito  c'lomc'iit  scheme's.  Finally,  the'  patch 
loads  dc't c'rminc'd  by  thc'sc'  intc'raction  models  can  hc'  applic'cl  to  highc'r  ordc'r  structural  modc'ls 
in  c'xactly  the'  same  mannc'r,  and  analogous  modc'ls  can  hc'  usc'd  for  multiple'  patch  jiairs  and 
more'  complc'x  gc'oinc't  ric's. 
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